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TRANSFER WITH UNIFORM NORMAL STRAIN 
by Robert G. Deissler 
Lewis Research Center 

SUMMARY 

Turbulence and longitudinal heat transfer for an axisymmetric accelerating or decel- 
erating flow were analyzed. The analysis was based on generalized two-point correlation 
equations, which were obtained from the incompressible Navier- Stokes, continuity, and 
energy equations. Viscosity, thermal conductivity, and the effects of mean strain were 
retained in the analysis, but the turbulence was assumed to be weak enough to neglect the 
triple correlation terms in the equations in comparison with the other terms. The cal- 
culated longitudinal turbulent heat transfer in the decelerating case was discussed in con- 
nection with the heat transfer at a stagnation point. 


INTRODUCTION 

This study is a continuation of work reported in references 1 and 2. In reference 1, 
weak locally homogeneous turbulence in a compressible flow with a mean longitudinal 
strain is analyzed. Reference 2 considers the turbulence in incompressible flow through 
a cone. The work in those references is, in turn, closely related to that of Pearson 
(ref. 3) and differs from earlier work mainly in the inclusion of the effects of viscosity 
and finite strain rates on the turbulence. 

In the present work, the effect of uniform mean velocity gradients dUj/dXp 
dt^/dXg, and dUg/dx^ on weak locally homogeneous turbulence and on longitudinal 
turbulent heat transfer is considered. Longitudinal turbulent heat transfer is of impor- 
tance, for instance, in considering heat transfer near a stagnation point. Locally homo- 
geneous herein means that the intensity of the turbulence does not vary appreciably over 
a correlation or mixing length. Shear stresses are assumed to be absent, and the flow 
is considered incompressible and axisymmetric. The turbulence portion of this work 
(no heat transfer) has been considered by Pearson (ref. 3). Additional results for the 


turbulence are given herein, inasmuch as Pearson gave results only for the accelerating 
case and did not include turbulent vorticities or spectra. Instead of a steady- state locally 
homogeneous but longitudinally varying turbulence, as considered herein, Pearson con- 
sidered a homogeneous time-varying turbulence. The two treatments, however, give the 
same results. 

The turbulence and turbulent heat transfer are assumed to be homogeneous in the 
transverse direction but only locally homogeneous in the longitudinal direction. The 
mean axial velocity is taken as uniform over a cross section. The turbulence is initially 
isotropic but is allowed to become anisotropic under the distorting influence of the mean 
flow. The basic equations required in the analysis will be considered in the next section. 


BASIC EQUATIONS 

The equations for locally homogeneous turbulent heat transfer will be considered 
first; the equations for the turbulence itself were obtained from reference 2. By writing 
the incompressible Navier-Stokes and energy equations at two points in the turbulent fluid, 
reference 4 (eq. (11)) shows that 
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Figure 1. - Vector configuration for two-point 
correlation equations. 


where the vector configuration for correlations between 
fluctuating quantities at points P and P’ is shown in 
figure 1. The quantity t is the fluctuating component 
of the temperature at P, u^ and uj are fluctuating 
velocity components at P and P', U k and Uj are 
mean velocity components, T is the mean tempera- 
ture at P, x^ and x^ are space coordinates, t is 
the time, p is the density, v is the kinematic viscosity, 
a is the thermal diffusivity, and p is the instantaneous 
pressure. The overbars designate correlations or aver- 
aged quantities. The subscripts can take on the values 1, 
2, or 3, and a repeated subscript in a term indicates a 
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summation. (All symbols are defined in the appendix.) From equation (17) in refer- 
ence 4, 
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For locally homogeneous turbulence and turbulent heat transfer, the turbulence is 
considered homogeneous over a correlation length, or the scale of the inhomogeniety is 
much greater than the scale of the turbulence. Thus, a quantity such as 

O p 

3 Tu!/a(x k ) a(x k ) in equation (3) will be negligible compared with 3 Tu!/9r, 9r,_. 

3 K m m J 

(A calculation for axially decaying turbulence without mean velocity gradients, ref. 5, 

fig. 3, implies that this is a good approximation except very close to the virtual origin 

of the turbulence. ) In general, for locally homogeneous turbulence, 3/ 3(x.) « 3/ 9r.. 

m 

Also, for that type of turbulence, the mean velocity and mean temperature will vary lin- 
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early over distances for which correlations are appreciable so that 3Ut/9x^ = 3Uj/3x k = 

3(IL) /3(xJ , UJ - U k = r 3(UJ /d(x ) , (XL + Uj/2 = (UJ , and 3T/3x k = 

1 m/ K m K K t K m / 4 m K K K m K 

3T /3(xJ . Finally, in order to make the set of equations determinate, the turbulence 


m 

is assumed to be weak enough to neglect terms containing triple correlations. The tur- 
bulence in a flow with large velocity or temperature gradients may not have to be as weak 
as that in a flow without mean gradients. The terms containing those gradients may be 
large compared with triple correlation terms, even if the turbulence is moderately strong. 
Equations (3) and (4) become, for steady state at a fixed point, 
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The case of uniform axisymmetric strain with no shear and with temperature gradient in 
the longitudinal direction is considered herein. Equation (5) for j = 1 and equation (6) 
then become 
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where ( l ) is not a tensor subscript. Equations (7) and (8) can be converted to spectral 
form by introducing the usual three-dimensional Fourier transforms defined as follows: 
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Then, using continuity and the inverse transform of equation (10) gives 


/•oo 

S ^^ elK; 


d/c 


(12) 


where iT is a wave number vector having the dimension 1/length and die = d«j di<^ dKn. 
Taking the Fourier transforms of equations (7) and (8) results in 
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where two equations have been combined into one by eliminating 

For axisymmetric strain, s ^ = s ( 3 )> an< l by continuity of the mean flow, 
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The turbulence is also assumed to be homogeneous in the transverse direction and changes 
only in the longitudinal or x^ direction, so that 
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To simplify the notation, let (U|) = U, (xj) = x, and S/^ = s in the remainder of the 
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paper. Then equation (13) becomes 
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where use was made of the continuity relation in the form x 2 y 2 + * 3^3 = Corre- 

sponding expression for <p^ and cp^ are given as equation (17) and (18) in reference 2. 
These expressions are 
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For uniform s = dU/dx, 
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Equations (16), (17), and (18) can be written in terms of c as 
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SOLUTIONS OF SPECTRAL EQUATIONS 


For solving equations (20) to (22), the turbulence is assumed to be initially isotropic 
(at c = 1), and 




(K 6.. - K.Kf 


(23) 


where Jq is a constant that depends on initial conditions (ref. 5, eq. (43)). For the 
initial condition on y^ (at c = 1), it is assumed that 

(Vi) = 0 (24) 

A 0 


That is, the turbulence-producing grid is assumed to be unheated, so that the temperature 
fluctuations are produced by the interaction of the mean longitudinal temperature gradient 
with the turbulence. Equation (23) appears to be the simplest condition that gives results 
that, at all values of x, reduce to those for isotropic turbulence as the mean strain goes 

to zero. The use of that initial condition for cp-- implies that Pearson's parameter 

2 ^ 
vKq/s approaches °°, where Kq is a characteristic initial wave number of the turbulence 

(ref. 3). Thus, the present results should be applicable for large kinematic viscosity, 
small initial turbulence scale, or small strain rate. The case VKq/s — 00 was not con- 
sidered by Pearson. 

Equations (20) to (22) can be solved by methods similar to those given in reference 6. 
Solutions of these equations subject to the initial conditions given in equations (23) and (24) 
are, in dimensionless form, 
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It can be seen that y|, and cp *. are functions only of c, k*, and Prandtl number. 

For Pr = 1, equation (25) can be integrated to give 
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In order to integrate over wave number space, spherical coordinates are introduced 
as follows: 


Kj = K COS 6 

(33a) 
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(33b) 

Kg = k sin cp sin 6 

(33c) 

For r = 0, equations (9) and (10) then become 

vi - r f'u dK 

(34) 

and 

tu7 = f r. d« 

j Jq J 

(35) 

where 

/•TT q 

\jy.. = / / cpiiK Sin 9 dcp d9 

1] Jo Jo 13 

(36) 

and 

y»7T y**27T rt 

1% = / / y.K sin 6 dcp dQ 

3 Jo Jo 3 

(37) 

The quantities and Tj are functions only of the magnitude of the wave number 

and represent spectrum functions that have been integrated over all directions in wave 

number space. 
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The expressions for the velocity variances u^, mm, u 2 , and Ug 

can be integrated 
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and c = U/Uq. 

Another quantity of importance is the turbulent vorticity tensor cckok. The vorticity 
spectrum tensor is given in reference 7 as 

n ij = ( 6 if 2 - K i K $<Pli ~ < 43 > 

As was the case for cp^, equation (43) can be integrated over all directions in wave 
number space to give 
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The vorticity tensor is then given by 


"‘"j ■ £ 


A ij 


dn 


(45) 


Calculated turbulent velocities, vorticities, temperature- velocity correlations, and 
spectra will be considered in the next section. 


RESULTS AND DISCUSSION 
Velocity and Vorticity Variances 

2 9 

Figure 2 shows turbulent velocity variances u^, Ug, and u^ plotted logarithmically 
in dimensionless form. Included in the plot is the curve obtained by solving equation (21) 
as though the effects of strain were absent by omitting the second to sixth terms. This 
solution gives 
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For an accelerating flow with uniform strain, the longitudinal component u^ de- 
creases more rapidly and the lateral components decrease less rapidly than they would 
if the effects of strain were absent. For large values of U/Uq, the lateral components 
reach a steady- state value as observed in Pearson's results (ref. 3). This result differs 
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from flow through a converging cone (ref. 2) where the increasing strain rate with dis- 
tance caused the lateral components to increase without limit as the apex of the cone was 
approached. The asymptotic equilibrium solution shown in figure 3 is given by 
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2 




(47) 


Thus, the solution represents a case in which the energy fed into the lateral components 
by straining action balances the energy dissipated by viscous forces. 

For decelerating flow near the virtual origin, both the longitudinal and the lateral 
components of the velocity fluctuations decrease in the direction of flow. For lower 
values of U/Uq, all components begin to increase as the effect of normal strain becomes 
greater than the effect of viscous dissipation. The region of increasing turbulent intensity 
in the decelerating case was not observed for a cone (ref. 2) where the strain rate s de- 
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Figure 


Longitudinal velocity ratio, c - U/Uq 

3. - Dimensionless variance of turbulent vorticity components for uniform incompressible strain. 


creased sufficiently with distance to allow the turbulence to decay. As U approaches 
zero in the present case, the turbulence components will tend to increase without limit. 
The assumption of local homogeniety will, however, tend to break down in that region, 
and the turbulence components will remain finite in a real situation. An increase in tur- 
bulent fluctuation in the decelerating flow near a stagnation point (in comparison with the 
free-stream fluctuation) has been observed experimentally in reference 8. 

The reasons for the trends observed in the turbulent velocity variances will become 

o 

clearer if the vorticity variances a>. plotted in figure 3 are considered. The dashed 
curves for no effects of strain were obtained from the equation 
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Consider first the accelerating case. Here the trends are opposite to those for the veloc- 
ity variances in figure 2; that is, the longitudinal vorticity component decays less rapidly, 
and the lateral components decay more rapidly, than they would if the effects of strain 
were absent. Thus, the turbulent vorticity tends to become alined in the flow direction. 
That alinement occurs first, because the longitudinal vortex filaments are strengthened 
by the stretching action of the mean flow, whereas the lateral filaments are shortened and 
thus weakened; and, second, because the mean strain rotates the vortex filaments which 
were originally oblique so that their axes tend to lie in the flow direction. The velocities 
associated with the turbulent vortex filaments will then tend to lie in the transverse direc- 
tions, in agreement with the curves for velocity variances in figure 2. As for the lateral 
components of the velocity variance, the longitudinal component of the vorticity variance 
approaches an equilibrium solution for large values of U/Uq in which the vorticity gen- 
erated by the mean strain balances that dissipated by viscous action. This solution is 
given by 



For decelerating flow at low values of U/Uq, the lateral components of the vorticity 
tend to increase, whereas the longitudinal component decreases more rapidly than it 
would if the effects of strain were absent. Thus, the vortex filaments tend to be alined in 
the transverse directions. This alinement occurs because the lateral vortex filaments 
are strengthened by stretching, while the longitudinal components are weakened since 
they are shortened, and because the axes of vortex filaments which were originally 
oblique are rotated toward the transverse directions by the stretching action of the mean 
strain in the transverse directions. With the turbulent vortex filaments mostly alined in 
the transverse directions, the velocities associated with them can be either in the longi- 
tudinal or the transverse directions. This explains why, for low values of U/Uq, the 
curves for both the longitudinal and the transverse components of the turbulent velocity 
variance in figure 2 increase in the flow direction, whereas in the curves for vorticity 
variance, only the lateral component can increase. The lateral stretching of the vortex 
filaments intensifies both the longitudinal and the transverse velocity fluctuations. 

Relative intensity ratios for turbulence components corrected to eliminate viscous 
decay are plotted in figure 4. For obtaining the ordinates in this figure, values of turbu- 
lent velocity variance with the effects of strain included (solid curves in fig. 2) are cor- 
rected to eliminate the effects of decay by dividing them by corresponding values for pure 
viscous decay (dashed curves in fig. 2). The result (after taking the square root) is di- 
vided by U/Uq to give intensity ratios that are relative to the local mean velocity. In 
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Longitudinal velocity ratio, c = U/Uq 


Figure 4. - Comparison of relative intensity ratios for turbulent components corrected 
to eliminate decay for several situations. 

addition to the present results for uniform strain in an incompressible flow, results for 
flow through a cone (ref. 2) and for uniform longitudinal strain in a compressible flow 
(ref. 1) are shown in the figure for comparison. The curves for uniform longitudinal 
strain in a compressible flow were obtained from equations (31) and (32) in reference 1 by 
noting that U g = Uq - sx Q and U g /y = Uq/Vq. The values for pure viscous decay were 
obtained by solving equation (23) in reference 1 with all but the first and last terms de- 
leted, and again using U g /i/ = Uq/Vq and U g = U Q - sXq. This solution gives (U/Uq) - " 7 ^ 4 

for the ordinate of the dot-dashed curve for i = 1 and (1 + U/Uq)^^(U/Uq) ^^ 4 /V2 for 
the ordinate of the dot- dashed curve for i = 2, 3. 

The curves for the lateral components (i = 2, 3) for accelerating flow are of partic- 
ular interest because, as shown in reference 2, the ordinates of those curves give 
approximately, for certain conditions, the Stanton number ratio St/StQ for the heat 
transfer between the fluid and a wall. In obtaining that relation, the normal strain is 
assumed to be so large that changes in the Stanton number along the flow are governed 
by normal strain rather than by shear. 

The curves for the lateral components (and thus for the Stanton number ratio) for 
accelerating flow in figure 4 indicate but a slight difference in the results for uniform 
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incompressible strain and for flow in a cone. That is, when plotted in this way, the re- 
sults at a given U/U Q for accelerating flow seem to depend but slightly on how the strain 
dU/dx varies along the flow. On the other hand the results for uniform longitudinal 
strain in a compressible flow lie significantly above the others. These results can evi- 
dently be explained by the fact that the stretching of the vortex filaments is more intense 
in that case, since the lateral compressive strain is absent because of the lack of a ra- 
dial flow. The dot-dashed curve might be related to heat transfer in a highly heated 
constant- area tube with fluid density changes along the length, whereas the other two 
curves are more closely related to nozzle heat transfer in which the effects of compress- 
ibility are small (ref. 2, fig. 8). 

Figure 5 shows the effect of uniform normal strain on dimensionless longitudinal 
turbulent heat transfer tu^*. Since ruj is divided by the temperature gradient b, the 
ordinates can be considered as representing the variation of longitudinal eddy conductivity 
with U/Uq. Results are given for Prandtl numbers of 0. 01 (liquid metals), 1 and 0. 7 
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(gases), and As Prandtl number decreases, the eddy conductivity decreases, appar- 
ently because a turbulent eddy in a high conductivity fluid, such as a liquid metal 
(Pr ~ 0. 01) gains or loses heat by conduction as it moves longitudinally and thus trans- 
fers heat with relatively low effectiveness (ref. 4, fig. 7). At a given Prandtl number, 
the trends in the curves are generally similar to those observed for the longitudinal com- 

2 

ponent of the velocity variance u^ in figure 2. A notable exception, however, is that 
tuj reaches a maximum at low values of U and then approaches zero at U = 0, where- 

as Uj became indefinitely large as U approached zero. The ratio of maximum to min- 
imum tuj is greater at low Prandtl numbers. 

The curves for decelerating flow in figure 5 illustrate the large increase that normal 
strain can produce in the longitudinal heat transfer between a body and a stream in the 
vicinity of a stagnation point when free-stream turbulence is present. The eddy conduc- 
tivity would decay to very low values (dashed curves) if the effects of strain were absent. 
(The time available for decay is quite large, since the fluid velocity becomes very small 
as the stagnation point is approached.) This increase in heat transfer is in agreement 
with the experiments in reference 9 and the analysis of reference 10. The increase is 
evidently produced by the lateral stretching of vortex filaments as assumed in refer- 
ence 10. Reference 10, however, considered only transverse vortexes, whereas the 
present analysis considers random vorticity in all directions. The present analysis does 
not, however, consider the damping effect of the wall at the stagnation point (effect of 
viscous diffusion), so that the increase in turbulent heat transfer due to normal strain is 
probably exaggerated here. In fact, an attempt to calculate the heat transfer near a 
stagnation point by assuming that the maximum in the curve in figure 5 (Pr = 0. 7 for 
gases) corresponds to the maximum eddy conductivity in the boundary layer, and that the 
minimum in the curve for i = 1 in figure 2 corresponds to the turbulence level in the 
undisturbed stream, gave (by dividing one ordinate by the other) increases in total heat 
transfer considerably higher than those observed experimentally. The results do indi- 
cate, however, that the combination of free-stream turbulence and normal strain (or 
lateral vortex stretching) can be an important factor in increasing the heat transfer in 
the vicinity of a stagnation point. 


Turbulence Spectra 

Dimensionless spectra of components of the velocity and vorticity variances are 
plotted in figures 6 to 8. The spectra show how contributions to the dimensionless mean 
velocity or vorticity fluctuations are distributed among dimensionless wave numbers. 
Plotted in this way, the curves for no strain (U/Uq = 1) reduce to a single curve that does 
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Longitudinal 


velocity ratio. 


Dimensionless wave number, k* - v" t (x - Xq)"^/ U jj 4, 
Figure 6. - Spectra of dimensionless longitudinal velocity variance u 


Longitudinal 
velocity ratio, 


Dimensionless wave number, /c* - v^lx - Xq^^kj Ujp 
Figure 7. - Spectra of dimensionless transverse velocity variance u 


Longitudinal 
velxity ratio. 


mm 


Dimensionless wave number, k* - v" fc (x - KjuQ 
Figure 8. - Spectra of dimensionless transverse vorticity variance w 






Figure 9. - Asymptotic equilibrium spectra for longitudinal ve- 
locity ratio U/Uq-oo. 



(b) Prandtl number, 0.01. 

Figure 10. - Spectra of dimensionless longitudinal eddy conductivity tu]*. 


not change with x. Thus, compari- 
son of the curves for various values 
of U/Uq shows how strain affects 
the spectrum at a given position in 
comparison with the spectrum at the 
same position with no strain. For 
instance, figure 8 shows that for de- 
celerating flow (U/Uq < 1) contribu- 
tions to the transverse vorticity 
occur at smaller wave numbers 
(larger vortexes) than they would for 
no acceleration. This trend seems 
to be congruous with the observation 
in the analysis of reference 10 that 
only the larger transverse vortexes 
are amplified by stretching in the 
neighborhood of a stagnation point. 
Figures 6 and 7 show that contribu- 
tions to components of the velocity 
variance also move to lower dimen- 
sionless wave numbers as velocity 
ratio decreases. For U/Uq > 1, the 
trends afe opposite to those for de- 
celerating flow, or contributions to 
the mean fluctuations move to higher 
dimensionless wave numbers. 

Figure 9 shows spectra corre- 
sponding to the asymptotic equilibrium 
solutions given by equations (47) and 
(49) for U/Uq — Note that the 
dimensionless quantities used in these 
spectra have been changed from those 
used in the preceding spectra in order 
to obtain finite dimensionless quan- 
tities for U — °°. The spectra in 
figure 9 are of interest because they 
show how contributions to the velocity 
and vorticity variances are distributed 
among wave numbers for a case in 
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which the energy fed into the turbulence by the mean strain exactly balances that dissi- 
pated by viscous action. Although the curves are in equilibrium at each wave number, 
there is not necessarily an equilibrium between production and dissipation at each wave 
number because energy can be transferred between wave numbers by the stretching of 
the vortex filaments by the mean velocity gradients as discussed in references 1, 2, 4, 
and 5. 

Spectra of the dimensionless longitudinal eddy conductivity tu^* are plotted in fig- 
ure 10 for Prandtl numbers of 0. 7 and 0. 01. The shifting of the curves to lower dimen- 
sionless wave numbers as U/Uq decreases is similar to that for the spectra of the longi- 
tudinal velocity fluctuations shown in figure 6. 

SUMMARY OF RESULTS 

For an incompressible accelerating flow with uniform strain, the longitudinal veloc- 
ity fluctuations decreased more rapidly, and the lateral fluctuations decreased less rap- 
idly in the flow direction than they would if the effects of normal strain were absent. 

For large values of velocity ratio, the lateral components were found to reach a steady- 
state equilibrium value, as observed in Pearson's results. This result differs from flow 
through a converging cone, where the increasing strain rate with distance caused the 
lateral components to increase without limit as the apex of the cone was approached. 

For decelerating flow at low values of velocity ratio both the longitudinal and trans- 
verse velocity fluctuations increased in the flow direction as the effect of normal strain 
becomes greater than the effect of viscous dissipation. A somewhat similar increase in 
velocity fluctuation in the decelerating flow near a stagnation point has been observed ex- 
perimentally. This region of increasing turbulent intensity in the decelerating case was 
not observed in an analysis for flow through a diverging cone, where the strain rate de- 
creased sufficiently with distance to allow the turbulence to decay. 

Many of the trends observed for the velocity fluctuation components could be ex- 
plained by the analytical result that the vorticity became alined in the longitudinal direc- 
tion for accelerating flow and in the transverse directions for decelerating flow. 

When the results for turbulent intensity were corrected for viscous decay and divided 
by local mean velocity, the transverse component for accelerating flow, which can be 
related to heat transfer between the fluid and a wall, was approximately the same for 
flow in a cone and for uniform incompressible strain. On the other hand the curve for 
uniform longitudinal compressible strain lay appreciably above the others, apparently 
because of the more intense vortex stretching for that case. 

The results for longitudinal eddy conductivity in a decelerating flow showed that 
normal strain could increase that quantity to values considerably above those which it 
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would have if the effects of strain were absent. Thus, free- stream turbulence with nor- 
mal strain (or lateral vortex stretching) could be an important factor in increasing the 
heat transfer in the vicinity of a stagnation point. 

Turbulent vorticity spectra showed that the turbulent vortexes in a decelerating flow 
tended to be larger at a given location than they would be for no deceleration. The eddies 
associated with the energy and with the longitudinal eddy conductivity also became com- 
paratively larger. Spectra were also obtained for the energy and the vorticity for the as- 
ymptotic equilibrium solutions at large velocity ratios. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 16, 1966, 

129-01-09-07-22. 
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APPENDIX - SYMBOLS 


b 

c 

C P 

H 

h 

J 0 

Pr 

P 

r 

St 

s, s (1) 
S (2 )> s (3) 


T, T_ 

’ m 

t 

Ui, (D^ 

Uj, (Uj)_, u 

U, 


m 


m 


longitudinal temperature gradient, dT/dx 

longitudinal velocity ratio, U/Uq 

specific heat at constant pressure 

defined by eq. (42) 

heat-transfer coefficient 

constant that depends on initial conditions 

Prandtl number, v/a 

pressure 

defined in fig. 1 

Stanton number, h/pc U 

r 

longitudinal velocity gradient, 3 Uj/9Xj = 3U/3x = 3(Uj) 



transverse velocity gradients, Sl^/Sxg = SOJg) /^(Xg) , 

/ / m/ m 

S(x 3 ) 

m 


mean temperature 
time 

mean velocity component 


mean longitudinal velocity component 


fluctuating velocity component 


U i U j 

x i> ( x i) 


m 


X, (xj 


m 


a 


defined by eq. (41), dimensionless velocity correlation 
defined in fig. 1 

longitudinal coordinate 


s(x - x Q ) 



thermal diffusivity 
defined by eq. (37) 
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'1 

r*i 

6 ii 

V 

e,cp 

K i 

< 

A ij 


A * 

A ij 


V 

p 

T 

TUi 


TU« 


<P 

<P 


xp*. 


n ii 




dimensionless spectrum function for longitudinal eddy conductivity, 

defined by eq. (10) 
defined by eq. (30) 

Kronecker delta 
defined by eq. (11) 
angular coordinates (see eq. (33)) 
wave number component 

defined by eq. (28) 

vorticity spectrum function, defined by eq. (44) 


v (x - x Q ) Tj 
J 0 U 0 


dimensionless vorticity spectrum function, A. 


i] 


■' V - * 0 ) 3 
U 0 J 0 


kinematic viscosity 
density 

temperature fluctuation 

longitudinal temperature-velocity correlation 


■' 6/ V - xnT- ™ 


3/2 


J u 3 / 2 b 

d o u o 

defined by eq. (9) 
defined by eq. (29) 

energy spectrum function, defined by eq. (34) 


xp.. 


■^(x - * 0 ) 2 
U 0 J 0 


defined by eq. (43) 
turbulent vorticity variance 
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w i w j 


w i w J 


, 7 /2 (x . Xq) 7/2 

uV2j 
U 0 J 0 


Subscripts: 

c corrected to eliminate viscous decay 

m at point P m (see fig. 1) 

w without effects of strain 

0 at virtual origin of turbulence where turbulent energy would be infinite (It is 

assumed that turbulence is isotropic at Xq and that strain begins to act there. ) 

1 in longitudinal or flow direction 

2, 3 in transverse direction 

Superscripts: 

’ at point P' (see fig. 1) 

* on dimensionless quantities 

- over averaged quantities 
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